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ABSTRACT 

We present an analysis of the spin-down parameters for 131 radio pulsars for which v 
has been well determined. These pulsars have characteristic ages ranging from 10 3 — 10 8 
yr and spin periods in the range 0.4-30 s; nearly equal numbers of pulsars have v > 
as v < 0. We find a strong correlation of v with i>, independent of the sign of v. 
We suggest that this trend can be accounted for by small, stochastic deviations in 
the spin-down torque that are directly proportional (in magnitude) to the spin-down 
torque. 
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1 INTRODUCTION 

Isolated pulsars exhibit fascinatingly rich timing behav- 
ior. The most dramatic variations are glitches, sudden 
increases in spin rate v, that are usually accompanied 
by increases in the magnitude of the spin-down rate 
(for a review of glitch properties see, e.g., Lyne, She- 
mar & Smith 2000, Krawczyk et al. 2003, and references 
therein). While large glitches have a well-defined signature, 
(Au, A;> ) = (+,—), microglitches exhibit all po ssible sig- 
natures llCordes. Downs fc Krause-Polstorr3ll98cfl . There is 
evidence that some isolated pulsars undergo "precession" 
(nutation) of their spin vectors (D'Alessandro & McCul- 
loch 1997; Stairs, Lyne & Shemar 2000; Shabanova et al. 
2001; Link & Epstein 2001; Akgiin, Link & Wasserman 
2006). Smaller, long-term variations in pulse phase, tim- 
ing noise, are seen in all pulsars. Timing noise has been 
largely attributed to high frequency random walks in the 
pulse phase, spin frequency or spin-down rate (see, for ex- 
ample, lloxntoneraLlIIizi ICordesI&Heifen^Iiiia) . Later 
analysis showed that the timing noise in some pulsars can- 
not be explained entirely in terms of high-frequency random 
walks, but that the spin behavior is due to discrete jumps in 
one or more of the spin parameters , possibly superimpose d 
on a high-frequency random walk llCordes fc Downslll985D . 
Some pulsars show quasi-periodic p hase variations over pe- 
riods of years (for examples see, e.g, iDowns fc Reichlvll983t 
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iD'Alessandro et aljll993T) . The term "timing noise" is per- 
haps a misnomer; to the extent that timing noise represents 
physical processes intrinsic to the neutron star, it is part of 
the signal. 

Upon construction of the pulse phase history of a pul- 
sar, 4>(t), higher order spin parameters (spin rate v and its 
derivatives v and v) are extracted from a Taylor series fit to 
the observed pulse phase, 



4>{t) = fa + v{t - t ) + -v(t - 



- t y 



+ - 6 Ht 



tof + 



(1) 



where fa is the phase at time to- Values of v determined in 
this way are typically orders of magnitude larger than the 
prediction of the vacuum dipole model and often differ in 
sign; such large values of v are thought to be highly noise- 
contaminated. [For example, Baykal et al. 1999, in a study of 
four pulsars that show significant quadratic trends in spin 
rate history, found instability in v consistent with a noise 
process]. In this connection, we note that determination of 
v is a difficult task since it is not a stationary quantity; val- 
ues deduced for a given pulsar generally depend on the origin 
of time for the fitting, the length of the fit and the partic- 
ular data set used (see, e.g., Hobbs et al. 2004, Fig. 7). As 
a result, there is considerable uncertainty in many measure- 
ments of v, and so most quoted values should be regarded 
as estimates. However, to the extent that v is determined 
by an underlying noise process that we would like to un- 
derstand, reliable estimates of v constitute a useful statistic 
with which to study the nature of timing noise in the pul- 
sar population as a whole. [We describe what we mean by 
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"reliable" below] . The idea behind the work we present here 
is to regard an estimate of v at a given time for a given 
pulsar as one outcome of the intrinsic "noise" process. We 
consider estimates of v in other pulsars as other realizations 
of a (possibly) similar underlying process. In this way we 
can quantify the strength of timing noise as a function of 
other spin variables: the spin rate, its derivative and the 
spin-down age. 

A complete characterization of timing noise has not 
yet been achieved, though trends of the strength of timing 
noise with various spin-down parameters have been identi- 
fied. Much previous work has used v as t he basic parameter 
with w hich to estimate the noise strength. ICordes fc Helfandl 
( 1980) quantified timing irregularities by defining an activity 
parameter: 



A = loa 



<J R {m,T) 
a R (m,T) C rc 



(2) 



where an(m,T) is the residual phase from a least-squares 
polynomial fit of order m over an interval of length T (all 
logarithms in this paper are base 10). To obtain (jfl(m, T), 
the measurement error to the rms residual is subtracted 
quadratically from the observed rms residual. Cordes & 
Helfand (1980) applied their definition of A for m = 2; in an 
analysis of 50 pulsars, they found evidence for correlations 
of A with spin period, period derivative, a nd spin-down age. 
More recently. I Arzoumanian et alJ lll994l) introduced a tim- 
ing noise parameter, characterizing the pulsar clock error 
caused by stochastic timing noise as 



A(/,=log( ±\i>\t 3 



(3) 



The statistic used bv I Arzoumanian et alJ l)l994h to quantify 
the strength of timing noise is Ag = A(t = 10 s s). The 
noise parameter is related to the activity parameter, A, by 
A = Ag + 0.42. In a study of 104 pulsars, Arzoumanian et al. 
(1994) found a correlation between As and p given roughly 
by 



A 8 = 6.6 + 0.6 log p 



(4) 



with a characteristic spread in As of a factor of ~ 3, inde- 
pendent of p. 

In this Letter we demonstrate the existence of a strong 
correlation between measured values of v and v. Nearly 
equal numbers of pulsars have v > as have v < 0. In- 
terestingly, we find that the correlations are independent of 
the sign of v. Because changes in the sign of v with time 
have been seen in some pulsars, we attribute the nearly 
equal numbers of pulsars with v > and v < as reflecting 
stochastic variations in v between two extremes over some 
time scale. In this interpretation, our results imply a corre- 
lation of the magnitude of timing noise (as measured by the 
characteristic magnitude of D) with the spin-down torque on 
the star (as measured by the observed i>). We suggest that 
the observed trend can be accounted for by small, stochastic 
deviations in the spin-down torque that are directly propor- 
tional (in magnitude) to the spin-down torque. 



2 AN ANALYSIS OF SPIN-DOWN 
PARAMETERS 

For this analysis, we have used the most reliable published 
values of the measured frequency second derivative v. These 
include the recentl y published ephemerides for 374 pulsars 
jHobbs et alj r2004) which have individual data spans of up 
to 34 yr and employ a new method to mitigate the effects 
of timing noise by whitening the timing residuals; the ro- 
tational parameters, particularly D, obtained in this way 
are significantly more precise than those obtained in earlier 
work in the sense that they give more stable timing solu- 
tions. While these v values can be used in a timing model 
to predict the pulse phase at an arbitrary time, Hobbs et 
al. (2004) warn that they are orders of magnitude larger 
than the values predicted by the dipole braking model and 
so may not reflect the physics of the pulsar braking mecha- 
nism (we offer an alternative point of view below). In gen- 
eral, the smaller \v\ are less stationary than larger values, 
and the relative errors can exceed 100%; we have therefore 
selected the subset of 127 pulsars for which the quoted er- 
rors in v are less than 10%. We also include four young 
pulsars whose v have been precisely measured, so that the 
Hobbs et al. (2004) de-w hitening process was not needed . 
These are 50531 + 21 llLvne. Pritchard fc Smith! Il988tl . 
5054 - 69 iCusumano Massaro fc Mined 120031). J1119 - 
6127 JCamilo et all2000D and 51509-58 feaspi et alll994l : 
lLivingstoT!^^t^lTl2005ir The data set shows a very large 
range in pulsar spin parameters. The characteristic ages are 
in the range ~ 10 3 — 10 s yr and the period 0.4-30 s. We 
have no millisecond pulsars in our sample. 

We consider correlations among the following spin pa- 
rameters: v, v, i>, the spin-down age r c = —vj2v and 
/>di P = 3u 2 /u. The definitions of r c and i>di P follow from 
the vacuum dipole model of secular spin evolution: v oc u 3 . 
Our results are plotted in Fig. 1, and the cross-correlation 
coefficients given in Table 1. The highest correlation (91%) 
is between |i>| and i>, with similarly strong correlations with 
!/dip and t c , and a weaker correlation with v. The trends of 
Fig. 1 can be summarized by the following fits: 



log |i>| = —11 + log 

log |£>] = -9.9 + 0.53 log i> dip , 

log |i>| = -8.9 - 1.2 log To- 



(5) 
(6) 
(7) 



This system of linear equations is self-consistent through 
the definitions of r c and i>di P - We have applied the analysis 
above to the full Hobbs et al. (2004) ephemerides of nearly 
400 pulsars, and we find the same trends found above for 
the sub-sample of 131. We note that: 

(i) The above trends are independent of the sign of v. 
The number of pulsars with v > is nearly equal to the 
number with v < 0. Both groups follow the same trends in 
magnitude of v (Fig. 1, second panel). 

(ii) \v\ is always greater than or equal to i>di P ; in many 
cases \v\ exceeds i)dip by factor of ~ 10 6 (Fig. 1, bottom 
panel). Only four pulsars fall on the line v = v&\ v : B0531+21, 
B0540-69, J1119-6127 and B1509-58. These pulsars appear 
on the upper right of Fig. Id. 

(iii) The correlation between \v\ and v shows generally 
less scatter than found by Arzoumanian et al. (1994) for As 
versus p. 
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Figure 1. The rotational and derived parameters versus v for 131 
pulsars. The (a) rotational frequency, (b) frequency derivative, (c) 
characteristic age (r c = —u/2u), and (d) the frequency second 
derivative predicted by the vacuum dipole model (i^ip = 3u 2 /u) 
Only four pulsars fall on the line v = u dip : B0531+21, B0540-69, 
J1119-6127 and B1509-58 (upper right of panel d). 



Aspect ii) implies very large braking indices, n — vi)/v 2 . The 
effective braking indices vary enormously among pulsars in 
the sample, between ~ — 1CP 6 and ~ 10 6 . This in no way 
implies a large spin-down torque (see below). We note that 
v is positive for the 12 pulsars with \v\~ti 10 -23 s -3 . While 
these pulsars comprise a small subset of our total sample, 
they do nearly follow the dipole prediction for v. For smaller 
values of \v\, there is a progression to values of v with both 
signs, accompanied by a marked deviation from the dipole 
prediction. These results are consistent with the notion that 
the values of v are noise-dominated for all but the largest 
values of \v\. Next, we show that >> £>di P might be due 
to small variations in the spin-down torque. 



3 DISCUSSION 

Theories of the physical origin of timing noise include star- 
quakes, internal torque fluctuations, external torque fluctu- 
ations and hybrid models that combine external torque fluc- 
tuations with magnetospheric changes that, in turn, affect 



Table 1. Cross correlation coefficients of pulsar rotational param- 
eters. We have not included cross-correlations among quantities 
that are defined in terms of each other, e.g., v and t c 
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the pul sar beam (for an exc ellent review of timing noise mod- 
els, see lD'Alessand"rclll997l) . We now discuss the possibility 
that the correlation between v and v is due to variations 
in the torque on the stellar crust. External torque fluctua- 
tions could arise from, for example, variations in currents in 
the pulsar magnetosphere (Cheng 1987a,b). Internal torque 
variations could result from vortex dynamics in the stellar 
crust (e.g., Alpar et al. 1986; Jones 1990). 

First, we point out that variability in the emission re- 
gion (through movement of the region above the star, for 
example) cannot, by itself, account for timing noise in pul- 
sars with relatively large values of \v\. A value of \v\ = 10 -24 
s~ 3 would lead to a phase variation of ~ 5 periods over 10 
years. For movement of the emission region to change the 
phase to this extent, the region would have to move many 
times around the star over time scales of years, without pro- 
ducing large changes in the pulse profile, an extremely un- 
likely scenario. In these pulsars, at least, it seems that the 
cumulative effects of timing noise over years must represent 
true variations in the rotational phase. [Short time scale 
variations, however, such as phase jitter, are likely due to 
processes in the emission region]. On the other hand, for 
pulsars with relatively small \i>\, we cannot rule out move- 
ment of the emission region as the dominant contributor 
to timing noise. A value of i> = 10 -28 s -3 would lead to 
phase variations of only < 10 -3 periods (0.4°) over 10 years. 
Since this phase difference is small compared to the typical 
beam width for these pulsars, the timing noise could be pro- 
duced by variations within the emission region. However, 
given that the correlation we have found between v and v 
holds over the nearly eight orders of magnitude in v that 
we have considered, a single process could be at work. We 
henceforth assume that timing noise represents variations 
in the spin rate of the star, and that these spin variations 
represent variations of the torque on the crust. We could be 
seeing variations in the spin-down torque on the star, or, in 
principle, variations in the internal torque exerted on the 
crust by the liquid interior. We consider the first process to 
be more promising. 

Let us write the total torque on the crust as 



Nt, 



N ad + N n 



(8) 



where JV s d represents the average spin-down torque and 
ATnoise is a stochastic (or, at least, largely non-deterministic) 
additional torque arising from fluctuations in the external 
torque, the internal torque, or a combination of the two. 
The observed v is 

V = !>sd + Z>noisc- (9) 

We expect that 1 2> no ise | is much smaller than | i> sc i | , while ob- 
servations require that |Aioise| is much larger than |i>di P j- 
This can happen if the characteristic time scale Tjv over 
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Figure 2. The torque parameter N = \u\Tn/\v\ versus spin-down 
rate. 



which the noise torque varies is short compared to the star's 
spin-down age. We can compare the characteristic noise 
torque to the spin-down torque with a dimensionless torque 
parameter, defined as: 



N 



(10) 



The torque parameter is shown in Fig. 2 versus where 
we take Tjv = 10 yr as a characteristic noise time scale for 
illustration (the exact value of Tjv is unimportant and, in 
any case, varies among pulsars). We point out two features 
from Fig. 2: 

(i) The torque ratio N shows no clear trend in v. Torque 



ratios of N <10~ 



account for the observed \u\(» u^ p ) 



in most pulsars. 

(ii) Because N is « 1, only small deviations from the 
spin-down torque are required to produce the observed val- 
ues of v. The magnitude of the noise torque is directly pro- 
portional to the external spin-down torque. [The scaling with 
the spin down torque can be seen in eq. JSJ]. 

The stellar magnetosphere is unlikely to be completely 
static, and so external torque variations are a natural ex- 
planation for our results. Whether or not internal torque 
variations applied to the crust by the liquid interior are suf- 
ficient to account for the range in v values presented here, 
is an interesting question. 

While completing this work, we became aware of a simi- 
lar analysis by Beskin et al. (2006), who find the same corre- 
lation of v with z> that we find. Our results are in agreement 
with theirs, but our interpretation is very different. Beskin 
et al. (2006) attribute the nearly equal numbers of pulsars 
with v > and D < to cyclic evolution of pulsar spin rates 
over time scales of hundreds of years. They suggest that 
their assumed spin variations could represent free-body pre- 
cession over very long periods, similar to that observed in 
PSR B1828-11 over a much shorter period (Stairs, Lyne & 
Shemar 2000). We emphasize, however, that the phase evo- 
lution over decades is essentially stochastic in many cases 
and quasi-periodic at best in other cases. Only a few pulsars 
show strongly periodic behavior. 
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